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We present cosmological constraints arising from the first measurement of the radial (line-of-sight) 
baryon acoustic oscillations (BAO) scale in the large scale structure traced by the galaxy distribution. 
Here we use these radial BAO measurements at z = 0.24 and z — 0.43 to derive new constraints 
on dark energy and its equation of state for a flat universe, without any other assumptions on 
the cosmological model: w — — 1.14 ± 0.39 (assumed constant), Q m — 0.24to'o5- H we drop the 
assumption of flatness and include previous cosmic microwave background and supernova data, we 
find w = -0.974 ± 0.058, Q, m = 0.271 ± 0.015, and tt k = -0.002 ± 0.006, in good agreement with a 
flat cold dark matter cosmology with a cosmological consant. To our knowledge, these are the most 
stringent constraints on these parameters to date under our stated assumptions. 

PACS numbers: 95.36.+X, 98.65.Dx, 98.80.Es 



Before recombination, the Universe was filled by a 
plasma of coupled photons and baryons. On horizon 
crossing, cosmological fluctuations produced sound waves 
in this plasma and when recombination occurred, 380 000 
years after the Big Bang, the distance covered by the 
sound wave was about 150 comoving Mpc. This is the 
so called sound horizon r s , also known as the BAO scale 
tbao- This signature can be found both in the cosmic 
microwave background (CMB) and large scale structure 
(LSS), so the baryon acoustic peak can be used as a stan- 
dard ruler in the Universe. The value of tbao depends 
on a few physical quantities, mostly the time to recombi- 
nation and the time of matter domination, both of which 
are well known from measurements of CMB temperature 
anisotropics, i.e. through parameters £l m h 2 and f^s/i 2 , 
also in good agreement with large scale structure and 
primordial abundance measurements. Thus the spec- 
trum of CMB fluctuations provide an accurate estimate 
of tbao = 153±1.9 Mpc, independent of the value of H Q , 
dark energy equation of state w or the curvature of the 
Universe p]. At the same time the CMB also provides an 
angular measurement of the BAO scale, which together 
with rsAO-, can be used for geometrical test, such as the 
measurement of curvature. 

A series of recent papers [H [31 [H [5] have presented the 
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clustering of Luminous Red Galaxies (LRG, see [6 ) in 
the latest spectroscopic Sloan Digital Sky Survey (SDSS) 
data releases, DR6 & DR7, which include over 75 000 
LRG galaxies and sample over 1 Gpc 3 /h 3 to z=0.47. The 
last paper [3] focuses on the study of the 2-point corre- 
lation function £(<r, tt), separated in perpendicular a and 
line-of-sight ir directions to find a significant detection of 
a peak at r ~ 110 Mpc/h (for H = 100 h Km/s), which 
shows as a circular ring in the a-ir plane. There is also 
a significant detection of the peak along the line-of-sight 
(radial) direction both in sub-samples at low, 2=0.15- 
0.30, and high redshifts, z = 0.40 - 0.47. The overall 
shape and location of the peak in the 2-point and 3- 
point function are consistent with its originating from the 
recombination-epoch baryon acoustic oscillations. This 
has been used to produce, for the first time, a direct mea- 
surement of the Hubble parameter H(z) as a function of 
redshift. This is based on a calibration to the BAO scale 
measured in the CMB [T]. The values of H(z) are then 
used in [5] to compare to other cosmological constraints 
that provide estimates for Hq. 

Here, instead of calibrating the BAO distance, we 
will use the direct dimensionless measurement as a red- 
shift scale (shown in Table |l| to provide cosmological 
constraints that are independent of Hq. This is some- 
how similar to what is done in Percival et al. [7] for 
the spherically- averaged (i.e. monopole) BAO detection. 
Both the 2-Degree-Field Galaxy Redshift Survey 2dF- 
GRS and SDSS spectroscopic redshift surveys have been 
used to constrain cosmological parameters at BAO dis- 
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tances U EH CCD EH ■ Photometric LRG catalogs have 
also been used to obtain cosmological constraints [L3lll4|. 

None of these earlier papers have constrained the 
baryon feature in the line-of-sight (LOS) direction. We 
will argue here that the radial BAO (hereafter RBAO) 
detection provides in practice fully independent con- 
straints from the monopole (or angular) BAO measure- 
ments. It is therefore important to check if such indepen- 
dent constrains agree with the flat concordance cosmolog- 
ical constant model (ACDM) and, if so, what are the new 
combined constraints when we add the new data to pre- 
vious results. Also note that in previous BAO detections 
the constraints come from either a fit to the overall shape 
of the correlation function or to the normalized wiggles in 
the power spectrum P(k), while in the LOS measurement 
°f r BAO m they only used the BAO peak location 
(identified by the position of maximum signal to noise 
around 90 — 120 Mpc/h) and not the shape. There are 
potential advantages and disadvantages of both methods. 
Using a fit to the overall shape is a more stringent test 
of the model, but puts more weight on the intermediate 
scales, where the statistical errors are smaller but where 
the modeling is more uncertain. Fitting the peak loca- 
tion is more model independent but, precisely because 
of this, provides little evidence for the physical origin of 
such peak. In the case of the analysis in [5], the phys- 
ical evidence comes from the agreement between model 
and data both on the radial direction (see Figs. 13-15 in 
[5]) and on the overall ^(tt, a) plane, which shows a clear 
BAO ring (with the right monopole contribution) that 
blends into the extended LOS peak (see Fig. 10 in [5]). 
A model with no radial BAO is disfavored at 3.2a. 

To work with fiducial comoving scales, [5] uses a refer- 
ence flat ACDM Hubble rate: 

H ref (z) = H oy /n m (l + z)* + 1 - Sl m (1) 

with Q m = 0.25 to determine fiducial values for tbao 
at two redshifts. Here, Azbao(z) are obtained from 
rBAo(z) in as Az BAO (z) = r BAO (z)H ref (z) /c. Note 
that the measurements of Azbao( z ) do not depend on 
Ho, since the dependences contained in tbao{ z ) an d 
H re f(z)/c cancel out. 

Previous BAO (e.g. [3 [8]) analyses have only looked 
at the monopole, which is the average of £(<r, tt) over 
orientations: 

6>(r)= / £(ff,7r)dM (2) 
Jo 

where r = \J o 2 + tt 2 and /i = n/r. The BAO peak 
found by [5] was located at radial tt ~ 110 Mpc/h for 
a < 5 Mpc/h, which corresponds to /i > 0.999. The con- 
tribution of radial modes with /x > 0.999 to the monopole 
at r = 110 in the above integral is less than 1% even if 
we take the amplitude of £(er, tt) to be 10 times larger in 
the radial tt direction than in the transverse a direction. 
Moreover, the covariance between radial and monopole 



correlations has been shown to be less that 10% on scales 
larger than 40 Mpc/h (see Fig. 2 of [15]). Thus, in prac- 
tice, we can regard here the BAO measurements from the 
monopole as independent from the RBAO measurements. 

Once the radial BAO scale has been measured [5] , there 
are several possible approaches to extract the cosmolog- 
ical parameters solely from this measurement. Percival 
et al. [7 j propose three possible ways. We have used two 
of them, that are described below. The third one uses 
the ratio of BAO scales at different redshifts. We have 
not used it, since we have the BAO scale at two different 
redshifts, so only one ratio can be constructed, and the 
degeneracy in the determination of cosmological param- 
eters is larger using this approach. 
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TABLE I: The BAO fiducial scale t B ao in the LOS direc- 
tion calculated with a flat reference H re f(z) cosmology of 
f2m = 0.25, for two redshift slices: z m is the respective pair- 
weighted mean redshift, and a st and a sy3 are the statistical 
and systematic errors on tbao (from [5]). Here we use the di- 
rect Azbao measurement, shown in the sixth column, which 
relates to the fiducial scale as Azbao = rBAoH re f(z)/c and 
is independent of the value chosen for H re f(z). The last two 
columns show the corresponding errors. 

We will derive cosmological constraints from the mea- 
sured values of Azbao shown in Table [I] which can be 
expressed as 



where r s is the sound horizon at recombination. It is im- 
portant to remark that the two measurements of Azbao 
are independent, due to the analysis strategy used in [5], 
which determines the BAO scale in two well separate in- 
tervals of redshift. 

Throughout this letter, we will assume a standard 
Friedmann-Lemaitre-Robertson- Walker cosmology hav- 
ing as free parameters the matter density f2 m , curvature 
fife, Hubble constant h (with H — 100 h km/s/Mpc), 
and the parameter w of the equation of state of dark 
energy, i.e. the ratio of its pressure to its density w — 
Pde I Pde- Unless otherwise specified, a constant w is 
assumed, and thus, the standard ACDM model is recov- 
ered for w = —1. We will start by assuming also a flat 
Universe, = 0. 

In order to constrain the dark energy parameters ^a(= 
1 — f2 m ) and w, we can take two approaches: 

a) We can express r s in (J3j as [2T] 

i 

c fi+^d da 

(4) 
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Data Set 
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+0.015 
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-0.0026 
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WMAP5 [J 


0.244 


+0.062 
-0.052 


-0.0049 


+0.0121 
-0.0061 


-1.13 


+0.37 
-0.39 


RBAO + 
WMAP5 + 
SNe [16] 


0.271 


+0.015 
-0.014 


-0.0021 


±0.0062 


-0.974 


+0.057 
-0.059 



TABLE II: Results from the cosmology fits with the radial 
BAO (RBAO) data and additional data sets assuming a con- 
stant equation of state parameter w. In the data set a) the 
BAO radial data are used, and the values and errors of fl m h 2 
and Q.f,h 2 from pQ are input into the computation of H(zi)r s . 
For the data set b) the BAO radial data is combined with the 
measurement of Ia(z*) from PQ. In all cases, the WMAP5 
measurements within the appropriate cosmology model (fiat, 
cosmological constant, etc.) are used, although the differences 
are minute. See text for details. The other rows combine the 
RBAO redshift scales presented in Table [I] with external data 
sets. All fits have values of x 2 /Ndof close to one. In particu- 
lar, the fit with RBAO, WMAP5 and SNe data within a non- 
flat constant-?/; cosmology (last row) has x 2 /N do f = 312/306. 
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FIG. 1: 68% CL constraints on the plane f2 m -u> when a flat 
Universe is assumed. In yellow, from the measurements of 
the RBAO scale plus flbh 2 and Q m h 2 . In pink, from RBAO 
plus l A and z* from WMAP5 pQ In blue, from the RBAO 
scale plus WMAP5 and in red from the RBAO scale, plus 
WMAP5 plus supernovae [16] . Results are consistent among 
them and compatible with a ACDM cosmology. 



b) Alternatively, we can use as external constraints 
the measurements by WMAP5 of the ratio Ia(z*) 
between the distance to the last scattering surface 
and r s (z*) and of z* itself, l A (z*) = 302.14 ± 0.87 
and z* = 1090.5 ± 1.0 with a ~ 40 % positive corre- 
lation, when assuming a flat CDM model with con- 
stant equation of state w. Again, Hq cancels out 
and we are left with only f2 m and w as unknowns, 
which we then proceed to determine. 



with 

_ 7 ft 2 (l + 0.2271iV e// ) 

aeq ~ n m h* [b) 

with Zd, the so-called drag redshift, written in terms 
of f^/i 2 and Q m h 2 (Eq. (3) in [T]), and take f^/i 2 = 
0.02273 ±0.0066 and n m h 2 = 0T329± 0.0064 from 
the five-year results of the Wilkinson Microwave 
Anisotropy Probe (WMAP5) [T] as external con- 
straints (we could have also used determinations 
from nucleosynthesis or LSS measurements), and 
fix fl 7 h 2 = 2.449 x 10~ 5 and N eff = 3.04. Note 
that the expression for AzBAo{ z i) m Eq. Q is in- 
dependent of Hq. Once the one-sigma WMAP5 
constraints on fl^h 2 and £l m h 2 arc included in the 
fit, the only remaining unknowns are ^2a(= 1 — ^m) 
and w. Since we have two independent AzsAo{ z i)i 
we will be able to determine them. 



In all cases, a frequentist \ 2 fit including correlations 
is applied to the data, using the publicly available Mi- 
nuit code |17j for minimization and contour calculation. 
Systematic errors in AzBAo{ z i) have been neglected, be- 
ing subdominant (see Table [I]) and relatively uncertain. 
Including them does not change the results significantly. 
The results obtained for Q m and w and their l-er er- 
rors, a(Q m ) and a(w), are shown in rows 1 and 2 of 
Table [ill Both determinations, a) and b), are consistent 
between them and in agreement with ACDM. The two- 
dimensional constraints at 68% CL can be seen in Fig. [l] 

Next, we have combined the results obtained from the 
radial BAO scale with the full CMB distance constraints, 
using the WMAP5 measurements of the shift parameter 
R(z*), the acoustic scale Ia(z*) and the redshift at decou- 
pling z* presented in Table 10 of pQ, and the covariancc 
matrix in Table 11 therein. We have also added the su- 
pernovae set compiled by Kowalski et al. TB], which can 
be found in [25]. We have taken their data covariance 
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FIG. 2: 68% CL constraints on the plane fi m -w when a 
non-flat constant-uj cosmology is assumed. In blue, from the 
RBAO scale plus WMAP5 [T] and in red from the RBAO 
scale, plus WMAP5 plus supernovae [16] . Results are consis- 
tent with a cosmological constant Universe with a precision 
around 5%. 

matrix without systematical errors. Adding systematics 
does not change the results qualitatively. The results of 
the combination assuming flatness can be seen in rows 
3 and 4 of Table [TT] We determine the equation of state 
parameter w, assumed constant, to be consistent with a 
cosmological constant to about 5%. 

If, instead, we assume the ACDM model and drop the 
assumption of flatness, we obtain the results in rows 5-8 
of Table |TT] which show consistency with a flat Universe 
at the 0.5% level. 

We can now keep the equation of state parameter 
w free and drop the assumption of flatness. The con- 
straints, shown in the last two rows of Table [TT] remain 
very robust when combining the three data sets. The 
two-dimensional Sl m -w contour at 68% CL is shown in 
Fig. [2] All results are in good agreement with ACDM. 
Moreover, the combination allows us to simultaneously 
determine w with a precision around 5%, and set the 
flatness of the Universe with a precision of about half 
of a percent. As a by-product, we obtain a Hubble pa- 
rameter h — 0.703 ± 0.020, consistent with most recent 



results in the literature. If instead of RBAO we use in 
the combination the BAO monopole determination in [8] , 
we find compatible results with errors that are 20-40% 
larger (see Table 6 in [IS]). On the other hand, if we 
combine RBAO and the monopole in [5], we find no im- 
provement in the error band with respect to the case with 
RBAO alone. 

Finally, we can drop the assumption of constant w 
and, instead, use the parametrization w(z) = wq + w a ■ 
z/(l + z) [T5J [IS]. Assuming a flat Universe, and com- 
bining RBAO with WMAP5 and supernova data, we get 
w = — 1.05±g:io and w a = 0.44±° with an 87% anti- 
correlation, consistent with ACDM (wq = — 1, w a = 0). 
We conclude that current data show no evidence for a 
time-dependent w, although the constraint on w a is still 
very weak. 

In summary, we have determined for the first time the 
cosmological parameters il m and w using the radial BAO 
scale, to be f2 m = 0.24±°;^ and w = -1.14±0.39. These 
results are perfectly consistent with the expectations of 
the ACDM cosmology. Moreover, when these results 
are combined with the cosmological distance constraints 
coming from the CMB [JJ and type-la supernovae |16| . 
we measure ft m = 0.271 ± 0.015, n k = -0.002 ± 0.006, 
and w = —0.974 ± 0.058. The result is in good agree- 
ment with a flat ACDM cosmology, and shows that the 
determination of the RBAO scale is consistent with all 
the other cosmological measurements. 

In the future, currently planned BAO surveys like 
Baryon Oscillations Sloan Survey (BOSS) [23] or Physics 
of the Accelerating Universe (PAU) [5D] will measure ra- 
dial BAOs with greater precision up to redshifts close 
to z = 1, producing more stringent constraints on the 
properties of dark energy. 
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